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We derive path-integral expressions for the second and third virial coefficients of monatomic quan- 
tum gases. Unlike previous work that considered only Boltzmann statistics, we include exchange 
effects (Bose-Einstein or Fermi-Dirac statistics). We use state-of-the-art pair and three-body poten- 
tials to calculate the third virial coefficient of 3 He and 4 He in the temperature range 2.6 — 24.5561 K. 
We obtain uncertainties smaller than those of the limited experimental data. Inclusion of exchange 
effects is necessary to obtain accurate results below about 7 K. 



I. INTRODUCTION 

Thermodynamic properties of fluids at very low temperatures are of significant interest. For example, the current 
International Temperature Scale [1] makes use of volumetric properties and vapor pressures of helium isotopes below 
the triple point of neon (24.5561 K); below the triple point of hydrogen (13.8033 K), the scale is based entirely on 
properties of 3 He and 4 He. The theoretical analysis of relevant properties at these conditions, such as the virial 
coefficients that describe the fluid's departure from ideal-gas behavior, is complicated by the presence of quantum 
effects. 

The inclusion of quantum effects in the calculation of virial coefficients was one of the first numerical applications 
of the Path-Integral Monte Carlo (PIMC) method. [2] In a series of pioneering works published in the 1960's, Fosdick 
and Jordan showed how to calculate the second and third virial coefficient of a monatomic gas using computer 
simulations. [3HS] Given the limited computational resources available at that time, they were able to calculate the 
third virial coefficient only in the case of two-body interactions, using a model potential of the Lennard- Jones form 
and assuming distinguishable particles (Boltzmann statistics). They argued that their method could be extended to 
include the proper quantum statistics, but they were able to compute exchange effects only in the case of the second 
virial coefficient. 

Recently, the exponential increase in computational power has enabled use of the path-integral method to calculate 
the properties of quantum degenerate systems, notably superfluid helium. [6 At the same time, progress in the 
computation of ab initio electronic properties of interacting atoms resulted in the availability of very precise two- 
and three-body interparticle potentials, at least for the lightest particles such as helium atoms [7HTT] or hydrogen 
molecules. [T2l U5] 

A natural application for these potentials is the calculation of virial coefficients. As is well known, the second 
virial coefficient depends only on the two-body potential, the third virial coefficient depends only on two-body and 
three-body interactions, etc. The second virial coefficient for a monatomic gas can be rigorously obtained at the fully 
quantum level from the calculation of the phase shifts due to the pair potential, and previous work has shown that 
a completely ab initio calculation of second virial coefficients for helium can have uncertainties comparable to and in 
many cases smaller than those of the most precise experiments. [T1HTB] 

In the case of the third virial coefficient, no closed- form solution of the quantum statistical mechanics problem is 
known. First-order semiclassical approaches have been derived [191120] and show that, in the case of helium, quantum 
diffraction effects result in significant modifications of the classical result, even at room temperature. However, there 
is no rigorous way to evaluate the accuracy or uncertainty of the semiclassical result, especially at low temperatures. 

In recent work, [3T] we extended the methodology pioneered by Fosdick and Jordan, deriving a set of formulae 
allowing a path-integral calculation of the third virial coefficient C(T) of monatomic species for arbitrary two- and 
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three-body potentials. Our results were limited to Boltzmann statistics (i.e., distinguishable particles) and we did not 
present results for temperatures lower than the triple point of neon (24.5561 K), which we deemed to be a reasonable 
lower bound so that exchange effects could be neglected. Nevertheless, we were able to compute the value of the third 
virial coefficient of 4 He with an uncertainty one order of magnitude smaller than that of the best experiments. 

Recent experimental results overlapping with our temperature range, |22l I23| although mostly consistent with our 
calculations, seemed to indicate a systematic deviation which the authors speculated could originate from our neglect 
of the proper quantum statistics of helium atoms. 

In this paper, we extend our computational methodology to calculate the quantum statistical contributions to the 
third virial coefficient, and compute C(T) for both isotopes of helium in the temperature range 2.6 — 24.5561 K, 
extending the temperature range considered in our previous work down into the range where exchange effects are 
important. We show that quantum statistical effects are significant only for temperatures smaller than about 7 K, 
and compare our results to low-temperature experimental data. 

In a subsequent publication |24j . we will present results covering the entire temperature range (improving on our 
previous results for 4 He at 24.5561 K and above) with rigorously derived uncertainties. We will also extend our 
methodology to include acoustic virial coefficients, and compare those calculations to available data. In the present 
work, our focus is on low temperatures and specifically on the effect of non-Boltzmann statistics. 



II. PATH-INTEGRAL CALCULATION OF THE VIRIAL COEFFICIENTS 

The second and third virial coefficients, B{T) and C(T) respectively, are given by [22] 

B(T) = -~{Z 2 -Z\) (1) 

C(T) = AB 2 (T) - ~ [Z 3 - iZ 2 Z x + 2Zf] , (2) 

where V is the integration volume (with the limit V oo taken at the end of the calculations), and the functions 
Zm are given by: 

Z 3 = A 9 f dld2d3 (mle-^^P^Im) (3) 

Z 2 = A 6 f dld2 (12|e-^ 2 ^7\ 2 |12) (4) 

Z x = A 3 / dl (l|e-^|l) = V, (5) 



where Hn is the iV-body Hamiltonian, ft = l/{k^,T), V is a permutation operator (multiplied by the sign of the 
permutation in the case of Fermi-Dirac statistics), the index ir^ runs over the 6 permutations of 3 objects (i.e., 123, 
132, 213, 321, 231 and 312), and ir 2 runs over the 2 permutations of 2 objects (i.e., 12 and 21). A = hj \f 7 2nvmk^F is 
the thermal de Broglie wavelength of a particle of mass m at temperature T. For the sake of conciseness, we denote by 
\i) an eigenvector of the position operator relative to particle i and by di (i — 1, 2, 3) the integration volume relative 
to the Cartesian coordinates of the i-th particle. Note that, in order to produce the molar units used by experimenters 
and in our subsequent comparisons with data, the right side of Eq. ([I ) and the second term in the right side of Eq. |2]) 
must be multiplied by Avogadro's number and its square, respectively. 

In the following, we will derive a path-integral expression for the calculation of the virial coefficients with Eqs. 
and ([2]). We perform the derivation in detail in the case of B(T) to establish the notation, and then extend the results 
to the more interesting case of C(T). 



A. Second virial coefficient 



In this paper, we adopt Cartesian coordinates to describe the atomic positions. This differs from the approach 
developed in Refs. 131 and |2"TI where Jacobi coordinates were used. This choice allows the exchange contribution to be 
computed in a much simpler manner than would be the case if Jacobi coordinates were used, especially in the case of 
three or more particles. 

From Eqs. |l]) and Q, it can be seen that there are two contributions to B{T). The first one comes from considering 
the identity permutation only, and takes into account only quantum diffraction effects. This is the only contribution 
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that gives a nonzero result at high temperatures, where the particles can be treated as distinguishable (Boltzmann 
statistics). 

The second contribution to B(T), which we will call exchange (xc), comes from the only other permutation involved 
in the definition of the quantity Z 2 above. 

The expression of these two contributions in Cartesian coordinates is: 



A6 

'zv 

A 6 



dridr 2 (rjT^I exp 



-P{K 2 + U 2 {\r 2 -r x \)) 



exp 



-0K 2 



\rir 2 ) 



B XC (T) = / drxdr 2 (r^lexp 



-0{K 2 + U 2 (\r 2 -r 1 \)) |r 2 n), 



(6) 
(7) 



where we denote by Km the total kinetic energy of N bodies and by U 2 (r) the two-body potential energy operator. 
The upper (lower) sign in Eq. (JtJ) corresponds to Bose-Einstein (Fermi-Dirac) statistics. 
Equations Q and ^ can be rewritten by using the Trotter identity 



e k 2+ u 2 = Um ( e K 2 /P e u 2 /p 

with a positive integer value of the Trotter index P. 

Following the procedure outlined in Ref. [2TJ one can then write i?Boitzmann(? 1 ) as 



-E>Boltz 



-2tt' 



r 2 dr (exp [-(3U 2 (r)] - l) 



where the two-body effective potential U 2 (r) is given by 



(8) 
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x 



1 



dr 
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+ U 2 (\r + Xl (r) -aj a (r)|)d.T 



where 



J ring 2 l 



p3/2 



A 3 



In the previous equations, we have defined Ax 



exp 



(10) 

(ii) 



(12) 



r^, , where r\!' is the coordinate of particle k (k = 1, 2) 



in the «-th "imaginary time slice" . These "slices" are obtained by inserting P completeness relations of the form 

1 = J drfdrf \r^rf){r^rf\ (13) 



between the factors e K2 ^ p and e U2 ^ p of the Trotter expansion of Eq. fl8J). We used the overall translation invariance 
of the system to remove the factor V in Eq. ^ and fix the r = slice of particle 2 at the origin of the coordinate 



system. We also denoted by x± and x 2 '' the coordinates of two ring polymers having one of their endpoints fixed 



.(<) 



at the origin (a4 1} = x\ L) = 0), and we introduced the variable r denoting the distance between the r = time slice 
of the two ring polymers. In the classical limit, where the paths x\(t) and x 2 (t) shrink to a point, the coordinate r 
reduces to the distance between the particles and one has U 2 (r) = U 2 (r). 



Note that the effect of the identity permutation is to set rY 



The path-integral formalism allows one to 



map the quantum statistical properties of a system with N distinguishable particles (Boltzmann statistics) onto the 
classical statistical properties of a system of N ring polymers, each having P beads (sometimes called imaginary-time 
slices), which are distributed according to the function F ring of Eq. (12 1. [55] The mapping is exact in the P —> 00 
limit, although convergence is usually reached with a finite (albeit large) value of P. 
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In the calculation of the second virial coefficient, Eq. ^ shows that the second virial coefficient at the level of 
Boltzmann statistics is obtained from an expression similar to that for the classical second virial coefficient, using an 
effective two-body potential. This effective potential, U 2 {r), is obtained by averaging the intermolecular potential 
U%{v) over the coordinates of two ring polymers, corresponding to the two interacting particles entering the definition 
of B(T). 

Equation ^ is equivalent to Eq. (19) of Ref. Ell The only difference is that the current approach uses Cartesian 
coordinates, and therefore we are left with an average over two ring polymers of mass m instead of one ring polymer 
of mass fi = m/2, corresponding to the relative coordinate of the two-particle system. The two approaches are of 
course equivalent, and in fact it can be shown that Eqs. ^ and (10 1 reduce to the form derived in Refs. 131 and |2"T1 
Equation (|9| is the same expression previously derived by Diep and Johnson for spherically symmetric potentials on 
the basis of heuristic arguments, [H] and later generalized by Schenter to the case of rigid bodies and applied to a 
model for water. [2"T] 

Equation (10) is actually the discretized version of a path integral, as shown in Eq. (11). The circled integral is 
defined as 



Dx exp 



1 r ph m 



hj 2 



dx(r) 



p-i 

lim / n dAxWF rto g(Ax< 1 >,...,A SB < p >) = l, 



(14) 



and it indicates that one has to consider all the cyclic paths with ending points at the origin, that is £c(0) = x(f3fi) = 0. 
The normalization of the path integral is also indicated in Eq. (14). 

We can perform on Eq. (JtJ) , describing the exchange contribution to the second virial coefficient, the same steps 
leading from Eq. (|6]) to Eq. fl9J) . The only difference is the presence of the permutation operator, whose main 

' +1 ^= To and ri P+1 " > — ri . In this case, defining X^ — r!'' 
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(15) 



(16) 
(17) 
(18) 



where we have defined A M = y/2A. The exchange contribution to the second virial coefficient is given simply as 
an average of the two-body potential taken on ring polymers corresponding to particles of mass /i = m/2. In the 
discretized version of the path integral, one has to consider 2P beads. In Eq. (17), we have used the overall translation 
invariance of the integral to remove the factor of V in the denominator. 

The effect of the various permutations can be visualized as generating paths with a larger number of beads, which 
are obtained by coalescing the ring polymers corresponding to the particles that are exchanged by the permutation 
operator. 
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B. Third virial coefficient 



Z 3 = A 9 / dld2d3 (123|e _/JHs 5^7^1123) 
Z 2 Z X = A 9 [ dld2d3 (123|e-^^ 2+f3 ^P 7r3 |123) 
Zf = A 9 [ dld2d3 (123|e-^ 3 |123), 



We now discuss the third virial coefficient, starting from the expression given in Eq. Since 4B 2 (T) can be 
calculated by the methods of the previous section, we concentrate on the second term, whose summands can be 
written as follows: 

(19) 
(20) 
(21) 

where Tj = — ^-V^ is the kinetic energy operator of particle i. 

We can simplify the expression in square brackets on the right-hand side of Eq. ^ by writing the three Z 2 Zi 
terms choosing each time a different particle for Z\ (in Eq. (20 1 we have chosen particle 3 as coming from Z\). After 
considering all the permutations of two and three particles, we end up with 6 + 3x2 + 1 = 13 terms building the 
term in square brackets of Eq. ([2]). It is useful to collect these 13 terms as follows: 

1. Term 1 (identity term): we sum together permutation 123 from Z 3 , the identity permutations from the three 
Z 2 Zi and the whole 2Z\ term. Adding 4£? Boltzmann (T), one obtains the Boltzmann expression for C(T), already 
discussed in Ref. [STJ In the present formulation based on Cartesian coordinates, the value C(T) in the case 
of Boltzmann statistics involves an average over three independent ring polymers, which correspond to the 
three particles. In the following, this contribution to C(T) will be referred to as CBoitzmann (T) and is made by 
1 + 3 + 1 = 5 of the 13 terms described above. 

2. Term 2 (odd term): we take permutations 132, 213 and 321 from Z 3 and the three exchange permutations from 
the Z 2 Z\ terms. These permutations are all odd, and we consider them with a positive sign (Bose-Einstein 
statistics). In the case of Fermi-Dirac statistics, this term has to be multiplied by an overall minus sign. All 
of these permutations correspond to configurations where two of the three particles are exchanged. The sum of 
these 6 terms will be referred to as C a dd{T). 

3. Term 3 (even term): we take the permutations 231 and 312 from Z 3 . These are the remaining two terms from 
the 13, and are both even permutations, hence the name. Both of these terms correspond to a cyclic exchange 
of the three particles, and their sum will be referred to as C even (T). 



Using these definitions, the full C(T), including quantum statistical effects, can be written as 

C(T) = CBoitzmann 

(T) ± C7 dd(T) + C cvon (T) + C B (T), 



(22) 



^B^+ABKT), 



(23) 



where the last term in the right-hand sum is given by 

Cb{T) = ±8£?Boitzi 

since the contribution of 4£?g oltzmann (T) to C(T) is already included in Csoitz 
upper (lower) sign corresponds to Bose-Einstein (Fermi-Dirac) statistics. 

Using the same procedure outlined above in the case of B(T), one can write the Boltzmann contribution to the 
third virial coefficient as 



(T). In Eqs. (22 1 and (23), the 



Boltzmann 
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(24) 
(25) 
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where -?v ing denotes the probability distribution of the path relative to particle k, as defined in Eq. (12 1. In Eq. (25 1, 

the three-body effective potential energy V 3 is obtained as an average performed over three independent ring polymers 
of the total three-body potential energy: 



V 3 (x,y,z) = U 3 (x,y,z) + U 2 (\x-y\)-\- 
U 2 (\x-z\) + U 2 (\y-z\), 



(27) 



where U 3 (x,y, z) is the non-additive three-body potential of three atoms. In Eq. (25) the total three-body potential 
energy for the Boltzmann contribution to the third virial coefficient is 



Vs(rur 2 ) = ^U 3 (r 1+ xf\r 2 + 4\x¥) 

i=l 

U 2 {\r l+ xf -r 2 -xf\) + 
U 2 (\r 2 + a-:, 



(0 
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.(0 



"2 ^ X 3^\)y 



(28) 



where the variables with superscript (i) denote the coordinates of three ring polymers with one of the beads at the 
origin. Notice that in passing from Eq. ^ to Eq. (24 1 we have used the translation invariance of the integrand to 
perform the integration over r 3 , which removed the factor of V in the denominator. As a consequence, the paths 
corresponding to particle 3 have their endpoints at the origin of the coordinate system (or, equivalently, the third 
particle is fixed at the origin when the classical limit is performed.) In the same limit, the variables T\ and r 2 
appearing in Eq. (26 1 reduce to the positions of particles 1 and 2, respectively, and one has V 3 (ri,r 2 ) = V 3 (ri,r 2 ). 

The term C dd(2j is obtained by exchanging the positions of two particles. This operation reduces the number of 
ring polymers to two: one having 2P beads, corresponding to the exchanged particles, and the other having P beads, 
corresponding to the remaining one. The odd contribution is given by 



C dd(T) 
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A 3 
2 3 / 2 
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-0H a 



exp 
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(29) 
(30) 
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<j) Vx exp 



0H 
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+ U 2 (\x(t + ph/2) - x(r)|)dr 



+ V 3 (x (r + ph/2) , x(t), r 3 + y(r)) dr 
1 , (31) 



where we have defined 



T7 odd 1 \ 
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(32) 
(33) 



The 2P variables X^ have been defined analogously to what has been done in Eq. ( 15 1. Notice that in the discretized 
version, the average defining the odd exchange term in Eq. (30) is performed over two different kinds of ring polymers: 
the first has 2P beads of mass m/2 and connects particles 1 and 2 whose coordinates are exchanged by the permutation 
operator, whereas the second - corresponding to the third particle of mass m - has P beads. 

A similar derivation holds for the even contribution to the third virial coefficient, which is given by 
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where we have defined 
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(*) y-(»+P) 



and 



yr(i+2P) _ r v«;_ j n |- ne discretized version, the even contribution to the 
third virial coefficient is an average over the coordinates of the 3P beads of a single ring polymer corresponding to a 
particle of mass to/3. 

Notice that, from a computational point of view, the evaluation of the exchange contributions to the third virial 
coefficient is much less demanding than the calculation of the Boltzmann part, which is given as an integral over 
the positions of two particles. In fact, the odd contribution is calculated as an integration over the position of one 
particle only, whereas the even contribution is given by a simple average over ideal-gas ring-polymer configurations. 
In particular, the full calculation of C(T) at the lowest temperature with 2.5 GHz processors required ~ 2400 CPU 
hours, only 15% of which was needed to calculate the exchange contributions. 



III. RESULTS AND DISCUSSION 
A. Details of the calculation 



We have calculated C(T) for both isotopes of helium with the path-integral method described above. We used the 
highly accurate two-body potential of Przybytek et al., |llj which includes the most significant corrections (adiabatic, 
relativistic, and quantum electrodynamics) to the Born-Oppenheimer result. We also used the three-body ab initio 
potential of Cencek et al., [10] which was derived at the Full Configuration Interaction level and has an uncertainty 
approximately one-fifth that of the three-body potential [5] used in our previous work. [2Tj 

We generated ring-polymer configurations using the interpolation formula of Levy. [3, 28] The number of beads 
was chosen as a function of the temperature T according to the formulae P = int[(1200 K)/T] + 7 for 4 He and 
P = int[(1800 K)/T] + 7 for 3 He, where int[x] indicates the integer closest to x. These values of P were enough to 
reach convergence in the path-integral results at all the temperatures considered in the present study. The spatial 
integrations were performed with the VEGAS algorithm [29] , as implemented in the GNU Scientific Library, [30] with 
1 million integration points and cutting off the interactions at 4 nm. The three-body interaction was pre-calculated on 
a three-dimensional grid and interpolated with cubic splines. The values of the virial coefficient and their statistical 
uncertainty were obtained by averaging over the results of 256 independent runs. 

First of all, we checked that our methodology was able to reproduce well-converged fully quantum B(T) calculations 
for helium, which were obtained using the same pair potential as the present work. [18] Our results agree within 
mutual uncertainties with these independent calculations, and confirm the observation, already made when analyzing 
theoretical B(T) calculations performed using Lennard- Jones potentials, that exchange effects are significant only for 
temperatures lower than about 7 K. |31j The exchange contribution to the second virial coefficient is negative in the 
case of Bose-Einstein statistics and positive in the case of Fermi-Dirac statistics, as one would expect. 



B. The third virial coefficient of He 



We report in Table [T] the values of the third virial coefficient of 4 He, together with the various contributions of 
Eq. (22), for temperatures in the range from 2.6 K to 24.5561 K, which is the lowest temperature studied in our 
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TABLE I: Values of the third virial coefficient of 4 He and its components at selected temperatures. The ± values reflect only 
the standard uncertainty of the Monte Carlo integration; see Ref . 1241 for complete uncertainty analysis. 



previous work. [21] The same data are plotted in Figure [l] where they are compared with the recent experimental 
measurements by Gaiser and collaborators. [221 123] 

More extensive comparison with available data over a wide range of temperatures will be presented elsewhere. |24j 
In Fig. [l] our results are plotted with expanded uncertainties with coverage factor k — 2 as derived in Ref. [M] the 
uncertainty at the same expanded level for the experimental results was estimated from a figure in Ref. 22 . 

First, we notice that exchange effects are completely negligible in the calculation of the third virial coefficient for 
temperatures larger than 7 K, where their contribution to the overall value is close to one thousandth of that of the 
Boltzmann part. This is analogous to what has already been observed for the second virial coefficient. 

When the temperature is lower than 7 K, the various exchange terms have contributions of similar magnitude and 
opposite sign, but their overall contribution to C(T) is positive at all the temperatures that have been investigated. 
The exchange contribution to C (T) is comparable to the statistical uncertainty of the calculation, which progressively 
increases as the temperature is lowered. 

In Fig. [I] it can be seen that our theoretical values of C(T) are compatible with those of recent experiments [22| [23] 
down to the temperature of 10 K. For lower temperatures, the experimental results are somewhat larger than the 
calculated values, even though agreement is found again for temperatures below 4 K, where C(T) passes through a 
maximum. 

C. The third virial coefficient of 3 He 

Similar behavior is observed in the case of the third virial coefficient for 3 He, whose calculated values are reported 
in Table |H| Also in this case the exchange contributions are of opposite signs, but their combined effect is to reduce 
the value obtained with Boltzmann statistics, which is the opposite trend to that observed for 4 He. 

The effects of the various contributions to the third virial coefficient, in both the Bose-Einstein and Fermi-Dirac 
case, are summarized in Fig. [2] for the representative temperature of T = 3 K. First, we notice that the largest 
contribution to the third virial coefficient comes from the Boltzmann term. The even exchange term has only a minor 
contribution, whereas the two remaining terms (C dd and Cb) have almost equal magnitudes and opposite signs. In 
the case of Bose-Einstein statistics, the contribution to C from C dd is negative, while that from Cb is positive; the 
opposite situation is observed in the case of Fermi-Dirac statistics. The overall sum of the exchange contributions is 
positive for 4 He and negative in the case of 3 He. 

The magnitude of each exchange contribution at a given temperature is significantly greater for 3 He; this reflects 
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FIG. 1: The third virial coefficient of 4 He. The black circles are the results of the present calculations, with error bars 
representing expanded uncertainties with coverage factor k — 2. The gray area shows the results of the recent low-temperature 
experiments by Gaiser and collaborators. [221 123] 



the larger de Broglie wavelength, which not only appears directly in the exchange terms but also affects the range of 
space sampled by the ring polymers. 

In the case of 3 He, the exchange contribution is significantly larger than the uncertainty of our calculations, at 
least at the lowest temperatures that we have investigated. Similarly to the case of 4 He, quantum statistical effects 
on C(T) contribute less than one part in a thousand for temperatures higher than 7 K. Even in the case of 3 He, we 
observe C(T) pass through a maximum, at a temperature around 3 K, which is 1 K lower than the temperature where 
C(T) reaches a maximum for the He isotope. 

There are only a few sources of experimental data for C(T) for 3 He. Keller [35] measured five pressure-volume 
isotherms at temperatures below 4 K; these were later reanalyzed by Roberts et al. [33 and meaningful values of C 
were obtained only for the two highest temperatures. A later analysis of the Keller data was performed by Steur 
(unpublished), whose equation for temperatures below 3.8 K was reported by Fellmuth and Schuster [33]. Some points 
were also extracted from volumetric data by Karnatsevich et al. [35] Recently, Gaiser and Fellmuth J3Hl [37] extracted 
virial coefficients from their measurements of two isotherms for 3 He with dielectric-constant gas thermometry. 

Figure [3] compares our calculated values to the available experimental data, where the error bars represent expanded 
uncertainties with coverage factor k — 2. Error bars are not drawn for our values above 5 K because they would be 
smaller than the size of the symbol. As was the case in our previous work, [5T] the uncertainty of our values of C (T) 
is determined by the statistical uncertainty of our Monte Carlo calculations (shown in Tables [i] and [n]) and by the 
uncertainty in the two- and three-body potentials. At the temperatures considered here, the statistical uncertainty is 
the dominant contribution to the overall uncertainty. The full uncertainty analysis is presented elsewhere. |24j 

For the experimental points, these expanded uncertainties were taken as reported in the original sources; we note 
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FIG. 2: The magnitude and sign of the various contributions to C(T) at T = 3 K. 



Temperature 




C 


CBoltzmann 


Codd 
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C B 




(K) 


(cm 6 


moP 2 ) 


(cm 6 


moP 2 ) 


(cm 6 mol 


- 2 ) 


(cm 6 mol -2 ) 


(cm 6 mol -2 ) 


2.6 


1338 


±29 


1857 


±28 


-1803 


±4 


-274.8 


±0.8 


-2047 


±2 


2.8 


1477 


±24 


1817 


±23 


-1164 


±3 


-167.9 


±0.6 


-1336.6 


±1.3 


3 


1480 


±17 


1712 


±17 


-760.2 


±2.2 


-105.7 


±0.4 


-886.4 


±0.9 


3.2 


1463 


±17 


1621 


±17 


-503.5 


±1.7 


-66.8 


±0.3 


-594.2 


±0.6 


3.5 


1395 


±13 


1487 


±13 


-277.2 


±1.2 


-34.89 


±0.15 


-333.7 


±0.4 


3.7 


1376 


±11 


1439 


±11 


-189.2 


±0.8 


-23.12 


±0.11 


-229.6 


±0.3 


4 


1303 


±9 


1342 


±9 


-107.0 


±0.5 


-12.69 


±0.07 


-133.43 


±0.16 


4.2 


1245 


±9 


1273 


±9 


-73.9 


±0.5 


-8.65 


±0.05 


-93.82 


±0.13 


4.5 


1173 


±7 


1190 


±7 


-43.7 


±0.3 


-4.86 


±0.03 


-55.94 


±0.09 


5 


1071 


±6 


1079 


±6 


-18.41 


±0.16 


-1.963 


±0.016 


-24.42 


±0.05 


6 


895 


±4 


897 


±4 


-3.56 


±0.06 


-0.353 


±0.005 


-5.143 


±0.013 


7 


772 


±3 


773 


±3 


-0.78 


±0.02 


-0.0784 


±0.002 


-1.196 


±0.005 


8.5 


645 


±2 


645 


±2 


-0.059 


±0.006 


-0.0087 


±0.0003 


-0.155 


±0.001 


10 


558.3 


±1.6 


558.3 


±1.6 














12 


475.5 


±1.1 


475.5 


±1.1 














13.8033 


426.2 


±0.8 


426.2 


±0.8 














15 


402.0 


±0.7 


402.0 


±0.7 














17 


369.6 


±0.5 


369.6 


±0.5 














18.689 


347.8 


±0.4 


347.8 


±0.4 














20 


333.4 


±0.4 


333.4 


±0.4 














24.5561 


297.8 


±0.3 


297.8 


±0.3 















TABLE II: Values of the third virial coefficient of 3 He and its components at selected temperatures. Note that the odd 
contribution Codd contributes with a negative sign to the overall value of the third virial coefficient C (see Eq. ( 22 1 ) . The ± 
values reflect only the standard uncertainty of the Monte Carlo integration; see Ref. [52] for complete uncertainty analysis. 
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FIG. 3: The third virial coefficient of ^He. 



that in some cases (notably Ref. [35} this appears to be merely the scatter of a fit and therefore underestimates the 
total uncertainty. 

Our results are qualitatively similar to the rather scattered experimental data. We are quantitatively consistent 
with the values based on analysis of the data of Keller, but values from the other experimental sources are more 
positive than our results. We note that a similar comparison for 4 He [24] . where the experimental data situation is 
much better, shows the C(T) values of Ref. [35] for 4 He to deviate in a very similar way not only from our results but 
from other experimental data we consider to be reliable. 



IV. CONCLUSIONS 



We used path-integral methods to derive an expression for the third virial coefficient of monatomic gases, including 
the effect of quantum statistics. We applied this formalism to the case of helium isotopes, using state-of-the-art two- 
and three-body potentials. 

We showed that exchange effects make no significant contribution to the third virial coefficient above a temperature 
of approximately 7 K for both the fermionic and bosonic isotope. This is the same behavior observed in the calculation 
of the second virial coefficient. For temperatures lower than 7 K, the sign of the contribution to C(T) from exchange 
effects depends on the bosonic or fermionic nature of the atom. In the case of 4 He, the exchange contribution to C(T) 
increases its value compared to the value obtained with Boltzmann statistics, although in our simulations the total 
exchange contribution has the same order of magnitude as the statistical uncertainty of the PIMC integration. In the 
case of 3 He, the exchange contribution is negative, and its magnitude is much larger than the statistical uncertainty. 

The range of temperatures that we have investigated covers the low-temperature maximum of C(T) for both 
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isotopes. The third virial coefficient of 4 He reaches its maximum close to 4 K, whereas in the case of 3 He the 
maximum is attained at a lower temperature. 

For both helium isotopes, the uncertainty in our calculated third virial coefficients is much smaller than that of 
the limited and sometimes inconsistent experimental data. For 4 He, we obtain good agreement with the most recent 
experimental results, except for some temperatures below 10 K. A full comparison with available experimental data 
for 4 He, including the higher temperatures of importance for metrology, will be presented elsewhere. [24] For 3 He, we 
are qualitatively consistent with the sparse and scattered experimental values; in this case especially our calculations 
provide results that are much less uncertain than experiment. In both cases, at the temperatures considered here, the 
uncertainty is dominated by the statistical uncertainty of the Monte Carlo integration, meaning that the uncertainty 
of C(T) could be reduced somewhat with greater expenditure of computer resources. 

We note two directions in which extension of the present work could be fruitful. One is the calculation of higher- 
order virial coefficients, which is a straightforward extension of the method presented here. This would be much more 
computationally demanding, but the fourth virial coefficient D(T) may be feasible, at least at higher temperatures 
where the number of beads in the ring polymers would not be large. Second, the method can be extended to calculate 
temperature derivatives such as dC/dT; such derivatives are of interest in interpreting acoustic measurements. Work 
on the evaluation of acoustic virial coefficients is in progress. |2"4"] 
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